A model of a dual-core matter-wave soliton laser 
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We propose a system which can generate a periodic array of solitary-wave pulses from a finite 
reservoir of coherent Bose-Einstein condensate (BEC). The system is built as a set of two parallel 
quasi-one-dimensional traps (the reservoir proper and a pulse-generating cavity) , which are linearly 
coupled by the tunneling of atoms. The scattering length is tuned to be negative and small in the 
absolute value in the cavity, and still smaller but positive in the reservoir. Additionally, a parabolic 
potential profile is created around the center of the cavity. Both edges of the reservoir and one edge 
of the cavity are impenetrable. Solitons are released through the other cavity's edge, which is semi- 
transparent. Two different regimes of the intrinsic operation of the laser are identified: circulations 
of a narrow wave-function pulse in the cavity, and oscillations of a broad standing pulse. The latter 
regime is stable, readily providing for the generation of an array containing up to 10, 000 permanent- 
shape pulses. The circulation regime provides for no more than 40 cycles, and then it transforms 
into the oscillation mode. The dependence of the dynamical regime on parameters of the system is 
investigated in detail. 

PACS numbers: 03.75.-b; 03.75.Lm; 05.45.Yv 

I. INTRODUCTION 

The concept of matter-wave lasers, which are intended to generate strong coherent beams of atoms, and many 
applications that such beams can find, are well known. It is commonly assumed that a Bose-Einstein condensate 
(BEC) may serve as a source of the matter-wave laser beam. Atomic lasing was considered in detail theoretically 
(see Refs. pj and references therein) and demonstrated in various experimental settings Q; the topic was recently 
reviewed in Ref. ||. Of special interest is a possibility to construct a laser operating in a pulsed (cmasi-soliton) 
regime, i.e., periodically releasing narrow localized pulses of coherent atom waves, as proposed in Refs. Q and [l| • In 
that connection, it is relevant to mention that (effectively) one-dimensional solitons in BEC with very weak attraction 
between atoms (in 7 Li) were created in well-known experiments 0. 

The first models of the pulsed matter- wave lasers assumed release of quasi-solitons from an elongated ("cigar- 
shaped") trap filled by a self-attractive BEC. This scheme, while probably realizable in the experiment, does not offer 
sufficiently effective means to control the operation regime, and does not secure generation of a very large number of 
pulses in a nearly periodic fashion. A more sophisticated scheme was very recently proposed in Ref. |7J , which assumes 
zero nonlinearity coefficient (i.e., zero scattering length of inter-atomic collisions) in a larger part of the elongated 
trap, and attraction between atoms (negative scattering length) in its smaller part, where the pulses are to be formed 
and released. The sign of the scattering length may be controlled and altered along the length of the trap by means 
of the Feshbach resonance (this mechanism was used for the experimental creation of the BEC solitons 6] ) . 

In this work, we aim to propose a different model of a pulsed matter-wave laser, in which the BEC reservoir and 
the pulse-generating cavity are separated. As shown in Fig. ^ below, they are elongated parallel traps with linear 
coupling between them, due to tunneling of atoms. Using the aforementioned possibility to control the scattering 
length by means of the Feshbach resonance, we assume that the interaction between atoms is (very weakly) repulsive 
in the reservoir, while in the cavity it is attractive (and weak too, see below). Additionally, the negative scattering 
length is assumed to be altered along the cavity's axis, which can be achieved using an appropriate configuration of 
the magnetic field responsible for the Feshbach resonance (we do not assume any time dependence of the scattering 
length or other parameters). 

Besides using the magnetic field, it was predicted theoretically and demonstrated experimentally 01 that the 
Feshbach resonance can be provided too by an adequately tuned optical field. Accordingly, a spatially modulated 
stationary distribution of the light intensity may also be employed to provide for necessary spatially inhomogeneous 
Feshbach-resonance configurations in the reservoir and cavity. 

It is relevant to mention that various soliton solutions in dual-core traps, with the linear coupling between the cores 
and opposite signs of the scattering length in them, were recently studied in detail |Tlj , following an earlier work which 
was dealing with a model of dual-core nonlinear optical fibers with opposite signs of the group-velocity dispersion 
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FIG. 1: A scheme of the proposed matter- wave laser composed of two nearly one-dimensional traps coupled by tunneling 
of atoms, with negative and positive scattering lengths in the cavity and reservoir, respectively. The opening at the right 
edge of the cavity (x = 1) depicts the valve releasing pulses into an outcoupling waveguide. Shading symbolically shows an 
instantaneous density distribution, |-i/;| 2 and |0| 2 , in both traps (larger density corresponds to more intense shading). The upper 
part shows the potential profile in the cavity, with the drop AU at x = 1 corresponding to the valve. 



(rather than nonlinearity) in the two cores |l2j |. 

In the model introduced below, both edges of the reservoir, and the left edge of the cavity are impenetrable to 
atoms, while the right edge of the cavity (a "valve" ) allows the release of pulses into an outcoupling atomic waveguide, 
see Fig. ^ The model was originally designed with an intention to provide for formation of a narrow quasi-soliton 
pulse in the cavity, that would periodically circulate in it, bouncing from the edges. Each time that the soliton hits 
the right edge (valve), a pulse is released into the outcoupling waveguide. After being slashed this way, the intrinsic 
solitary pulse is supposed to replenish itself in the course of the subsequent cycle of the circulation by collecting atoms 
tunneling from the reservoir. Accordingly, the initial condition is taken with a large number of atoms in the reservoir, 
and a small number in the cavity. 

As described in detail below, numerical simulations demonstrate that the circulation regime outlined above may 
be observed, but it is not a really stable one: the intrinsic pulse would perform no more than 40 circulations, quickly 
coming to a halt and getting broad. Nevertheless, a very robust regime of the periodic release of pulses can be found 
in a large parameter region. In this regime, a broad "lump" stays immobile in the cavity, performing periodic cycles of 
stretching and compression. When its right wing reaches the valve, in the course of each cycle of the vibrations, a new 
outcoupling pulse is released. The respective loss in the number of atoms in the cavity is compensated by the influx of 
atoms tunneling from the reservoir. The pulse-release cycles repeat very many times, till the reservoir gets essentially 
depleted. We stress that the model is a passive one, i.e., it supports the stable pulse generation regime without any 
externally applied active control, such as periodically opening and shutting the valve by some clock signal. 

The paper is organized as follows. In Section 2, the model is formulated. Basic results, which demonstrate 
both unstable circulation-based and stable vibration-based regimes of the periodic generation of solitary pulses, are 
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presented in Section 3. In Section 4, we analyze the dependence of the operation of the matter-wave laser on its 
parameters, and give an estimate for the actual number of atoms in outcoupled pulses generated by the laser. The 
paper is concluded by Section 4. 

II. THE MODEL 

According to what was said above, the model of the matter-wave laser is based on a system of linearly coupled Gross- 
Pitaevskii equations (GPEs) for the wave functions tp(x,t) and </>(x,i) of atoms in the (effectively) one-dimensional 
cavity and reservoir, with opposite signs in front of the nonlinear terms. We scale h, the atomic mass m, and the 
nonlinearity coefficient at x = in the cavity (in physical units, the latter is Airh 2 a/m, where a is the s-wave scattering 
length of atomic collisions around x = 0) to be 1. Then, the coupled equations are cast in a normalized form, 

iipt = —^xx - 0- + g%) 

+/(s-L/2)V-«6 (1) 
i( f>t = --(f>xx + e\(j>\ 2 (/) - Kip, (2) 

where the coordinate x takes values < x < L (i.e., L is the length of both the cavity and reservoir). In accordance 
with what was said above, we assume a small positive scattering length (weak repulsion between atoms) in the 
reservoir, e > 0, while in the cavity the interaction is attractive, with the strength linearly increasing with x (g > is 
assumed), following a corresponding profile of the static magnetic or optical field that controls the scattering length 
in the cavity via the Feshbach resonance. 

The coefficient g will be small (see below), but it is needed to create an effective potential slope in the cavity: as the 
intra-cavity pulse gets "heavier" due to the agglomeration of atoms tunneling from the reservoir, the term proportional 
to g creates a gradually increasing effective potential which pushes atoms to the right (since g is positive). In the 
aforementioned circulation regime, this effective (nonlinear) potential pushes the intra-cavity localized pulse to the 
right. After hitting the right edge (at x — L), where the "valve" is placed (see details below), the pulse loses a 
part of its mass, generating the outcoupled soliton. Therefore, the force induced by the nonlinear potential drops, 
allowing the intra-cavity pulse to move in the opposite direction after the bounce from the right edge, until it will 
bounce (without loss) from the left edge, x = 0. This way, the nonlinear term proportional to g may support periodic 
circulations of the internal pulse (although, as said above, in the domain of the model's parameter space explored in 
this work the circulations are not really stable, eventually switching into the regime of vibrations; the latter regime is 
possible with g = as well). 

Further, we assume that an external parabolic potential acts in the cavity, 

U po t(x) = /(x-L/2) 2 , (3) 

with the center set at the middle of the cavity. It is well known that such a potential corresponds to a magnetic or 
optical trap applied to the condensate . It will be shown below that the results are not sensitive to the exact shape 
of the intra-cavity potential; in fact, essentially the same stable vibration regime can be obtained in the model with 
no potential. Note that Eq. J2J for the reservoir does not include any potential, as it is not necessary to generate the 
dynamical regime sought for. We also explored a model including a potential in the reservoir, but it did not affect 
the operation regime in any remarkable manner. 

The linear coupling between the cavity and reservoir is accounted for by the positive coefficient k, which determines 
the tunneling time for atoms. Keeping e,g,f and k as free parameters of the model, the remaining scaling 
invariance of Eqs. and J2Jl makes it possible to set L = 1, which is fixed below. 

As said above, the model assumes that the reservoir at its both edges, and the cavity at its left edge are closed by 
impenetrable lids, which gives rise to the respective boundary conditions (b.c), 

(j)(x = 0) = (p(x = 1) = tp(x = 0) = 0. (4) 

The pulses (solitons) are to be released into an external waveguide from the cavity at its right edge, x = 1 (recall we 
have set L = 1). For this purpose, a valve is set at x = 1, assuming that the corresponding b.c. for the -0 field is 
linear. Actually, the single possible form of such a linear b.c. is 



ip x (x = 1) = iqip(x = 1), 



(•5) 
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with a positive constant q (its physical meaning is discussed below). Indeed, using the continuity equation for the 
total density, p = \i[>\ 2 + \<fr\ 2 , and current, j = (i/2) (ip*4> - ip x i>*) + (*/ 2 ) {4>* x <t> ~ <Px<P*) = ji + hi m tne coupled 
GPEs © and ©, 

| = >^* (6) 

integrating Eq. © over the length of the system, < x < 1, and taking into regard that ji(a; = 0) = ]2{x = 
0) = J2(x = 1) = pursuant to b.c. Q), one can derive a depletion equation for the solution's norm N, which is 
proportional to the total number of atoms in the system: 

N = [ (\4>(x)\ 2 + \0(x)\ 2 )dx = N 1 +N 2l (7) 
Jo 

~ = jt(x = l,t) = q\ij(x=l,t)\ 2 = E(t) > (8) 

where b.c. 10) was substituted for ip x in the expression for ji(x = 1). The latter result is what might be expected: the 
instantaneous rate at which the density is released through the valve is proportional to the local density \ip(x = 1)| , 

which justifies the adoption of b.c. ©. Initial values of the norms n[ ^ and of the ip and <p fields, and the b.c. 
constant q are supplementary control parameters of the model, in addition to the above-mentioned set of (e, <?, /, k). 
Below, it will be demonstrated that really important parameters are e, n, and q. 

The physical meaning of b.c. © can be readily understood. Indeed, it implies that the local gradient of the phase 
of the wave function ip at the point x = 1, which is proportional to the atom's momentum, is fixed to be q. This 
can be realized by assuming that, to the right of the point x = (in the outcoupling waveguide), the cavity potential 
drops by AU — q 2 /2 (recall we have set m = K= 1). The drop AU must be essentially larger than the potential and 
kinetic energy of atoms to the left of the point x — 1 (in the cavity) . The latter condition, which amounts to 

U(x - 1) = //4 « q 2 /2 (9) 

[see Eq. JjJJ], provides for the possibility to neglect small deviations in the kinetic energy of the released atoms from 
q 2 /2. We notice that a soliton-releasing valve in the matter- wave laser model recently proposed in Ref. includes a 
similar element (an effective potential step). 

The form of the function E(t), as defined in Eq. ©, determines the shape of the density pattern (actually, an 
array of pulses) released into the outcoupling waveguide, the same way as the temporal shape of the input signal at 
the point z = determines the shape of a temporal soliton or solitonic array propagating along the coordinate z in 
the optical fiber ^4|. Note that, as follows from Eq. 10), the limit cases of q = and q = oo correspond, respectively, 
to a reflecting mirror and impenetrable lid, i.e., ip x {x = 1) = and ijj(x = 1) = 0. In either case, the system does not 
release anything. In addition to yielding the final result in the form of E(t), Eq. ©, with dN/dt computed directly 
from numerical data, provides for a means to monitor the accuracy of simulations. 

Finally, we note that the model does not include an explicit form of the outcoupling waveguide (the extension of 
the cavity beyond x = 1), which is responsible for the eventual shaping the released pulses into solitons. Actually, the 
soliton- formation problem in a uniform waveguide has been studied in detail before (see, e.g., Ref. |5| and references 
therein). 



III. OPERATION REGIMES OF THE PULSED MATTER- WAVE LASER 



In the numerical analysis of the model, Eqs. Q an d © were solved by means of a finite-element pseudospectral 
method, which implies division of the integration domain into several finite elements. Solution in each element is 
approximated by a polynomial, and a pseudo-spectral method, based on Chebyshev collocation points M 

, is used in 

each element. We have checked that a numerical algorithm based on domain stretching, rather than domain division 
into finite elements 0], produces the same results. 

In cases when the stationary version of the GPEs, i.e., a system of two linearly coupled ordinary differential 
equations, had to be solved (see below), the finite-element pseudospectral method reduces them to a set of algebraic 
equations, which were treated by means of the Newton's method. In the general case, when Eqs. and (J2J) are 
partial differential equations, the time derivatives in the equations were approximated by the finite difference as per 
the Crank-Nicholson scheme. This, together with the finite-element pseudospectral method, yields a set of nonlinear 
algebraic equations, that were also solved by the Newton's algorithm. 
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FIG. 2: Typical examples of the initial real profiles of 'iz(x) (a) and (b) for g\ = 0.02, e = 0.1, and k = 0.5, used in 

simulations of the coupled equations Q and @. In all the cases, the norm of the 0-component is fixed (in the normalized 
units) to be N2 = 363.4, while the norm iVi of the ^-component varies as shown in the figure. 



In all the cases considered, the integration domain < x < 1 was divided into four equal elements, with a fifth-order 
polynomial used in each one. The procedure leads to a set of 48 algebraic equations, and it was checked that increase 
of the polynomial's order did not change the results in any tangible aspect. We used the time step of At = 0.001, 
checking that smaller time steps gave virtually the same results (while larger At would generally be insufficient to 
produce accurate results). 

Simulations started with initial conditions which, by themselves, were stationary solutions to Eqs. Q and (J5J, in 
the form of {ip(x,t), 4>(x,t)} = exp (— ifit) {&(x), &(x)} with real fy(x) and <I>(a;). These stationary solutions were 
found, in a numerical form, for given values of g, e and k. While doing so, b.c. |J5J at x = 1 was temporarily replaced 
by ip = [as said above, this formally corresponds to q = oo in Eq. JSJ], since a stationary solution is not possible 
otherwise, due to the nonzero current ji across the point x = 1, see Eq. (|HJ). Further, both equations Q and used 
to generate the stationary solutions were modified by including ad hoc potentials linear in x, as the so generated initial 
configurations were found to quickly produce stable operation regimes, even though the actual potentials in Eqs. 
and 10) are different [recall that Eq. (J2J) has no potential at all]. A set of typical initial configurations obtained this 
way is displayed in Fig. |21 

Starting from such initial configurations, simulations produced two types of dynamical regimes. The first one 
features a relatively narrow pulse which performs shuttle motion, circulating in the cavity. Each collision of the 
pulse with the valve at the right edge of the cavity gives rise to a quasi-soliton released by the matter-wave laser. 
However, as mentioned above, the circulation regime was never found to be stable (although we cannot claim that 
the investigation of the system's parameter space was exhaustive, as it is difficult to perform complete exploration of 
the seven-paramater space) . The pulse would perform no more than 40 circulations (typically, fewer) , and this regime 
would transform itself into a stable vibration regime. 

An example of the transition from circulations to vibrations is displayed in Fig. [3J Actually, the transition may be 
much shorter at other values of the parameters. In most cases, as shown below, the vibration regime sets in directly 
from the initial configuration, without going through the transient stage of circulations. 

The instability of the circulation regime can be understood. Indeed, one may consider a small fluctuation which 
makes the norm of the pulse released into the outcoupling waveguide slightly smaller than in the established regime, 
hence the norm remaining in the cavity soliton becomes slightly larger. Consequently, the strength of the effective 
nonlinear potential, that tries to push the soliton to the right, drops (as a results of the bounce from the valve) 
by an amount which is a bit smaller than in the unperturbed state, and the soliton will therefore slide to the left 
slower than in the absence of the perturbation. As a result, spending more time in the cavity, the soliton will absorb 
more atoms in the course of the next cycle of the circulations, and will return to the right edge with a still larger 
norm. The mechanism by which a fluctuation slightly increasing the soliton's norm leads to its still larger increase 
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FIG. 3: Transition of the circulation regime into a vibration one, shown in terms of a three-dimensional plot for the density 
\ip{x, t)\ 2 in the cavity (a), contour plots of the density (b), and the outcoupling rate E(t) (c). The parameters are g = 0.02, 
e = 0.1, / = 3, K = 0.5, and q = 50. The initial (normalized) numbers of atoms in the reservoir and cavity are iV 2 (0) = 360 and 
JV< 0) = 12. 

(or, conversely, a fluctuation decreasing the norm causes its further decrease) implies the instability. It may happen 
that the circulation regime can be stabilized in a model including an active element (i.e., forced circulations might be 
stable), but in this work we focus on the more fundamental passive model. 

In the regime of vibrations, a broad pulse stays at the center of the cavity and performs a very large number of 
internal oscillations. At a stage of the vibration cycle when the pulse spreads out, its right wing hits the valve and 
generates an outcoupling localized matter- wave packet. A typical example of such a stable regime is displayed in Figs. 
Q] and O The second panel in Fig. 0] shows (for a long interval of time) the temporal dependence of the outcoupling 
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FIG. 4: A typical example of the establishment of a stable vibration regime in the cavity, shown in terms of the outcoupling 
rate E(t), which determines the shape of the pulse array generated by the matter- wave laser [E is defined as per Eq. JSJ]. 
(a) The initial stage; (b) beginning of the established regime; (c) a fragment from panel (b), showing 15 cycles in detail; (d) 
a single typical cycle (shown in terms of obviously rescaled variables t* and E*), together with a simplest analytical fit to it. 
In panels (a) and (b), contour plots in the lower subpanels (ii) additionally illustrate the spatiotemporal evolution of the atom 
density, \i/j(x,t)\ 2 . Parameters are g = 0.02, e = 0.10, / = 3, k = 0.5, and q = 200. The initial normalized amounts of matter 
in the reservoir and cavity are Nn = 363.4 and Ni = 6. In the established regime, the period is T = 0.135, with the duty cycle 
of 50%. 



rate E(t), as defined by Eq. ©. As said above, this dependence actually determines the shape of the pulse array 
released by the matter- wave laser. 

Figure O depicts vibrations of the density distribution in the cavity, for the same case as displayed in Fig. As seen 
from these figures, after a transient period, a very robust regime (with some residual long-period modulations) sets 
in. With the maximum value of the release rate E in this regime E max ~ 0.4, the oscillation period T — 0.135 and the 
duty cycle (the share of the period within which the value of the density at the right edge of the cavity, \i/i(x = 1, t)\ 2 , 
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FIG. 5: The evolution of the density distribution in the cavity, \ip(x, t)\ 2 , for the same case as shown in Fig. 0] 

exceeds half of its maximum value) being 0.5 imply that the number of cycles which are expected before depletion of 
the reservoir will become appreciable (say, N will drop to two thirds of = 363.4) can be easily estimated: 



(2/3)ivf 
(T/2) 



10, 000 cycles. 



(10) 



Direct verification of this prediction requires too long simulations (up to t ~ 1500). An estimate for the actual number 
of atoms in the generated pulses (which may range between 10 and 1000) is given in the next section. 

It is relevant to stress that, in the cases presented in Figs. |21and0][51 with / = 3 and q = 200, the condition 
which substantiates the use of b.c. ©, is satisfied with a huge margin. In fact, this condition holds equally well in all 
the cases where stable operation of the matter-wave laser model was observed. It is also relevant to mention that, in 
physical units, the value of q ~ 100 corresponds, for 7 Li, to the velocity ~ 1 cm/s of the released atoms, the respective 
kinetic energy being ~ 0.1 nK, on the temperature scale. 

It is relevant to compare the depletion rate in the stable operation mode [as given by Eq. l|5)l] with the rate at which 
the cavity and reservoir exchange the matter. A generic example of this comparison is displayed in Fig. As is seen, 
the exchange rate is much higher than the speed of depletion. It is interesting to note that the long-period beatings 
in the dependence of Ni(t), which are observed in El practically do not manifest themselves in the oscillations of the 
outcoupling rate E(t) [nor in the evolution of Ni + N2, as the beatings in Ni(t) are almost exactly compensated by 
anti-phase beatings in A^(t)]. In other words, the beatings do not affect the quality of the pulse array generated by 
the matter- wave laser. 

The line showing Ni + N2 vs. t in Fig. seems to have finite thickness due to fast oscillations of E(t), such as 
those in Fig. 0] Its comparison with the respective dependence predicted by averaging of the depletion equation © 
(a bold dashed line) is also shown in Fig. |SJ An error due to the averaging may account for a small discrepancy with 
the actual evolution of iVi + N2. 



IV. DEPENDENCE OF THE OPERATION REGIME ON PARAMETERS 



The matter- wave laser model proposed above depends on several parameters. It is necessary to identify those which 
are really important to the stability of the operation mode, and which are not. First of all, simulations demonstrate 
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FIG. 6: The evolution of the norm (amount of matter) in the reservoir, Ni (the rapidly oscillating dashed curve), and the 
total norm, Ni + N% (the continuous finite-thickness line) , in a stable operation regime. The bold dashed curve displays the 
prediction following from the depletion equation (JHJ, averaged over rapid oscillations of E(t). Parameters are the same as in 
Fig. H 

that there is virtually no dependence on the potential's strength / in Eq. (Q: at least within the interval of 1 < / < 3, 
all dynamical characteristics of the matter-wave laser model display very little variation, if other parameters are fixed 
(taking typical values of physical parameters given below, we conclude that, in this parameter region, the range of 
1 < / < 3 actually corresponds to the range of the longitudinal trapping frequencies 0.1 — 1 Hz). The dependence 

on the initial amount of matter in the cavity, N[ , is very weak too, at least in the interval of 5 < n[ ^ < 20. Also 
negligible is sensitivity of the established regime to the parameter g in Eq. QJ. In fact, additional simulations have 
demonstrated that virtually the same stable oscillation regime sets in if we take Eqs. and @ with / = g = 0. 

A more important role is played by the constant q in Eq. © (recall it determines the momentum of atoms in the 
outcoupling pulse), and the linear-tunneling constant n that determines the strength of the linear coupling between 
the reservoir and cavity in Eqs. Q and J5J. Figures [7] and [S] illustrate effects of the variation of these parameters on 
the operation regime. This is shown via the change of the time dependences of the outcoupling rate E(t) and total 
amount of matter remaining in the system, Ai + N^. 

The essential dependence of the dynamical regimes on the parameters q and k is quite natural. In addition to this, 
and somewhat unexpectedly, the system demonstrates strong sensitivity to the value of the nonlinearity coefficient 
e in Eq. which is proportional to the positive atomic collision length in the reservoir. In fact, stable operation 
regimes, providing for the generation of very long periodic arrays of the outcoupling pulses, are possible in a relatively 
narrow interval of values of e, if the other parameters are fixed. The dependence of the dynamical states on e is 
illustrated, by displaying segments of the established dependence E(t), in Fig. |5J Panel (a) in this figure shows that, 
at e = 0.08, the generated array of pulses has a high amplitude, but with shallow dips between pulses, so that they 
are not well separated. The same panel demonstrates that the optimum shape of the array of well-separated pulses 
is achieved between e = 0.09 and 0.12 (which implies the positive scattering length in the reservoir being ~ 10 times 
smaller than the absolute value of the negative scattering length in the cavity). The situation in the latter interval 
is additionally illustrated by a set of curves in panel (b) of Fig. [5] Physically, the necessary fine-tuning of e in the 
reservoir may be achieved by means of the Feshbach resonance. 

Undoing normalizations that cast the coupled GPEs in the rescaled form of Eqs. JTJ and (J2J, it is easy to derive a 
relation between the actual number of atoms, A p h ys , and the norm Ai of the tp function in the normalized units: 



(11) 
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FIG. 7: The influence of the variation of the boundary-condition parameter q (the momentum of atoms in pulses released by 
the matter- wave laser), see Eq. on the time dependence of the outcoupling rate (a), and on the depletion dynamics (b). 
Except for q, other parameters are the same as in Fig. [I] 




where r is the transverse radius of the cigar-shaped trap, which plays the role of the cavity, a and L being, respectively 
(as defined above), the Feshbach resonance-controlled atomic scattering length in the cavity and the cavity's length, 
both taken in physical units. For the cases displayed in Fig. El a characteristic value of the norm of an individual 
released pulse, in the normalized units, is Ni ~ 0.1. Assuming the aspect ratio of the cylindrical trap to be Lj (2r) ~ 10, 
and the actual length L ~ 100 fxm, Eq. (|TT)l yields the number of atoms per outcoupling pulse between 10 and 1000, 
if the scattering length is kept at a low level, between a ~ —0.1 nm and a ~ —0.01 nm (a very small value of \a\ 
is necessary to prevent collapse in the cavity trap). These values may be sufficient for physical applications. In 
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FIG. 9: The effect of variation of the nonlinearity coefficient e in the reservoir on the operation regime of the matter-wave laser 
model. Other parameters are the same as in Fig. 3. 

particular, the threshold number of 7 Li atoms, necessary for the creation of a stable three-dimensional solitons in an 
optical lattice, was recently shown to be 60 |17|. The number of atoms ~ 100 may also be sufficient for direct optical 
detection of the matter-wave pulse. As for the total number of atoms originally stockpiled in the reservoir, the above 
estimates, including Eq. (|10l) . yield the range of iVo ~ 10 5 — 10 7 for it. 

V. CONCLUSION 

The objective of this work was to propose a model of a matter-wave laser which can provide for stable periodic 
generation of a sequence of solitary-wave pulses from a reservoir containing a large amount of matter in the form of 
a coherent Bose-Einstein condensate. The system includes two parallel cigar-shaped traps, viz., the reservoir and the 
work cavity, which are coupled by the tunneling of atoms between them. The scattering length of atomic collisions 
is tuned to be positive and negative in the reservoir and cavity, respectively. Both ends of the reservoir, and the left 
edge of the cavity are closed by lids. Solitary pulses are released through a "valve" at the right edge of the cavity, 
which is described by the linear boundary condition (J5J , and may be implemented as a potential step separating the 
cavity and the outcoupling waveguide. Two different regimes of the operation of the matter-wave laser were identified: 
circulations of a narrow solitary pulse in the cavity, and vibrations of a broad standing lump. Only the latter mode 
is stable, while the circulation regime spontaneously rearranges itself into the vibration mode. Dependence of the 
stability and characteristics of the pulse-generating regime on parameters of the model was explored. The regime is 
sensitive to the boundary-condition parameter q and strength k of the linear coupling between the core and reservoir, 
and especially sensitive to the value of the nonlinearity coefficient (i.e., scattering length) in the reservoir. 

The bottom line is that the model can guarantee stable generation of a pattern consisting of up to 10 4 permanent- 
shape pulses, each containing between 10 and 1000 atoms, under typical experimental conditions. It is relevant to 
stress that the parameter space of the model, which is huge (seven-dimensional, as the full set of the parameters 

includes e, g, /, k, g, and the initial norms N^), is far from being fully explored, so still more promising regimes may 
be hidden in the model. 
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